3. Image Enhancement Using Filtering 
In the Spatial Domain 


* Chapter 3 (Image Enhancement in the 
Spatial Domain) 


Filtering for: 
- Removing noise 
- Sharpening image 
- Detecting edges 


1-D Linear Time Invariant (LTT) 
Systems 


1-D signal x(n) passes through a system T: 
y(n) = T [x(n)] 
A sysyeait jaa (τὴ) =a, y (1) + a; y, (3) 


A system igimeymyarignt: 


An LTI system cape f'aapletely 
characterized by its, impulse response: 


1-D convolution = 2 210117- κ) = PÅ I(k)x(n- k) 


2-D Linear Shift Invariant (LSI) 
oystems 


2-D c 


2-D signa n. m) passes through a system 
TÉ un,n)=T 222 ر‎ n») 


A زا۳۳۵9‎ PPE (qn) -ay(,n)* a y (1,2) 
A syst RAS- Salis Inyariapt; - m,n- m) 


The LSI system can be completely 


(point sp MES DUE HOUSE response 
Ham) = δ Y xk, k)n- kn- k) 


o je 


=x*h 
kd, gl * Ma 


mm, 
oe a A 

E m ^ 
καὶ - c 

= اسح سم‎ Tu 
ΞΞ ' τας πῶ 
z ge S a 
ee 
e 


= 
Ο 
D 
= 
© 
> 
= 
© 
© 


A Typical 
2-D 


h(n1,n2)[3 Δ 
hà 


m 


A 2-D 

Convoluti 

on 

Example ie sone ^ | 


Ky 


(f) 


ο Ed 
4 10 10 6 
4 10 10 6 
Y 3 35 2 


A Separable 2-D System 


Given an M x M system ry , it is 
PERB) SG (3) where 
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Row-wise 
convolution 


Convolution 
with 2-D Separable 
Filter 
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A Separable 2-D System 


(1): for a ΗχβΗ , we perform 1-2 convolktigh of 
ho». and 


(2): for a fixed  , we perform 1-D convéintibn of 
ha). and 


Total # of multiplications: (for N > > M, 
reduction facto —M/2) 


~NM 2 -—2N'M 


Some Other Useful Properties 


A 2-D LSI system n D) is stable, if 
b S n.n) > مه‎ 
n=” n=” 

A 2-D LSI syste n, n,) is causal, if 


An, n) = n «0 or n, «0 
Typical Lowpass Filter: done by 
averaging. 


Typical Highpass Filter: done by 
differencing. 


2-D Convolution in Image 
Processing 


* Most filters used in image processing 
have 

coefficients symmetrical with respect to 
the 

center. In these cases, in performing the 
2-D 

convolution, we don't need to flip h(n,, 
n»). 


e When the mask is partly outside the 
image, 
we often assume the values of the pixels 
outside the image have the same values 


1-D Filters 


Lowpass Highpass 


Frequency domain 


û 0 0 
(a) (b) (c) 


Spatial domain 
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A 2-D Lowpass F 
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Image after Lowpass Filtering 


Lowpass Filtering for Noise 
omoothing 


Most image energy is in its low frequency 
components (high spatial correlation among 
neighboring pixels), while most image noise is 
wideband (e.g., white Gaussian noise) -- calls for 
lowpass filtering (LPF) (it cannot remove the 
noise in the low frequency range though !). 


LPFs cause blurring, which is critical for 
images with sharp edges, but not critical for 
images of smooth contrast. 

Another use: remove blocky effect created in 
lossy image compression. 
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Some Typical LPFs 


Averaging LPFs: 
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Lowpass 
Filtered 
Images 


E Average 
filter of size: 
3X9; DXO, 
7x7, L95xl5; 
and 25x25. 


Median Filtering 


To avoid blurring in noise removing. Median 
filtering is a nonlinear process especially 
useful for reducing impulse, or salt-and- 
pepper noise. 


It preserves edges. 


Median filtering: with a sliding window 
passing through the image, choose the 
median (middle) intensity value -- nonlinear. 


Sorting can be done recursively. 
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Median Filters Preserve Edges 
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Averaging Mask vs. Median 
Filter 


abc 


FIGURE 3.37 (a) X-ray image of circuit board corrupted by salt-and-pepper noise. (b) Noise reduction with a 
3 X 3 averaging mask. (c) Noise reduction with a 3 X 3 median filter. (Original image courtesy of Mr. Joseph 
E. Pascente, Lixi, Inc.) 
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A Highpass Filter 
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Highpass Filters 


= A basic highpass spatial filtering; 4 
-1 8 -1 


-1 -1 -1 
E Note that the sum of coefficients is Ὁ! The DC term of 


H(w1,w2) is zero, i.e., the global contrast of the images is 
reduced. 


* Results need level-shift and scaling or clipping to remain in 
[0,1 -1 |. 
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Derivative Operators 


An averaging is analogous to "integration", 
"differentiation" can be expected to have 
the opposite effect: sharpening. 


Differentiation (gradienf(ooP? 
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Derivative Operators 


Typically, we want to average over several 
RENDRE o i.e., 
2 αἰ füy*1n- D- fa-12- 1) 


on 
+ fín*1n)- fn- 1,1) 
+| f(n +1,n +1(- f(n-1n +1) 


These result in Prewitt operators: 


-1 01 -1 -1 -1 
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Prewitt Operator 


-10 1 
-101 operate on 
-10 1 
۶0 - ل - وبا‎  f(um-D f(u*in-1 
f(n- 1n) f(n,n) f(n *1n) 


f(n - 1,1 +1 f(n, n, * 1) {η *1,n, 11) 


=| f(n+1n,-1- f(n-in-1) 
+ f(n+1n)- fn- 1n) 
+| fn +12 +1)- fn- 12 +1) 
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Derivative Operators 


How about emphasize the center row more: 


EN ~| f(a +1, - Ώ- fia-1n-1) 


+ 2 f(n+1n)- f(n- 1,1) 
+| f(n*«1,-41)- f(n- 1 +1) 


These result in Sobel operators: 
-1 O 1 -1 -2 -1 
-2 ο 2 0 0 0 
-1 0 1 1 2 1 
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Diagonal Derivative Operators 


We can also deal with both derivatives 


Εμ) TUR) n). 


F ~ f(n,n)- f(n + با‎ n +1) 


This results in Robert operators: 


1 0 0 1 
0 -1 -1 0 
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Effects of 1st & 2nd Derivatives 


First and Second 


a b 
C 
FIGURE 3.38 
(a) A simple 
image. (b) 1-D 
horizontal grav- 
level profile along 
the center of the 
image and 
including the 
isolated noise 
point. 
(c) Simplified 
profile (the points 
are joined by 
dashed lines to 
simplify 
interpretation). 
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Second Derivative —1 0 0 00 1 0 6-126 0 0 


Laplacian Operators: 2nd 
Derivative 
"f(n, "fp, 
VA. n rip ,O Fay ry 
on on 
Difference Eq. approximation of Laplacian: 


-- ~ f(n*1)- f(n,n) 


«fin *1n)- f(1,n))- | fü,n)- fü-1n)| 
= f(n*1,n )- 2f(n,n (+ fQ- LB) 

V fin,n) = fn *1,n)* f(1 - 1,1) + 

f(n,n *1)* f(n,1,- D- Af(n,n) 


δ΄ f(n, n) 
η΄ 
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Laplacian Operators 


0 1 0 1 1 1 -1 2 
1-41) 1-81 2 - 4 
0 1 0 1 1 1 2 2 


The edges are indicated by zero crossings 
(which offer more reliable edge location), 
instead of large gradients. 


Enhancement Using Laplacian 
Operators 


f(x y)- Ν΄ f(x گر(‎ centebaplaciarpeffisnegati 


gx y) -| f(x, y) + ۱۷ f(x yif centebaplaciarpeffis positiv 


g(x y = f(x y)- [f(x+1, y+ f(x- 1, (+ f(x y*1) f(x y- 1) 
- Af(x y]-5f(x pD-Lf(x*1 (+ f(x- 1, * f(x y+1+ f(x y- 1] 
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peg 

FIGURE 3.40 

(a) Image of the 
North Pole of the 
moon. 

(b) Laplacian- 
filtered image. 
(c) Laplacian 
image scaled for 
display purposes 
(d) Image 
enhanced by 
using Eq. (3.7-5). 
(Original image 
courtesy of 


NASA.) 
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High-Boost Filtering 
| -1 -1 -1"1-1 -1 -1[ [0 0 0 
= With A21, -1 A+8 -1-2|-1 8 -1 +0 AO 


a b 

c d 

FIGURE 3.43 

(a) Same as 

Fig. 3.41(c). but 
darker. 

(a) Laplacian of 
(a) computed with 
the mask in 

Fig. 3.42(b) using 
ΑΞ 0. 

(c) Laplacian 
enhanced image 
using the mask in 
Fig. 3.42(b) with 
A = 1.(d) Same 
as (c), but using 
A = 1.7. 
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High-Boost dew ue 


vq gs «οὐ 
E With A21, 


A-1) (original 
. Highpass 


(a) 
original 
(b) 
A=1.1 
(c) 


A=1.15 


(d) 
A=1.2 


Image after High-Boost 
Filtering 


